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Abstract: The lepton flavor violation /x 67, r — > 67, r ^ /X7, /x ^ 3e, r 3e, 
T 3fi and fi — e conversion in Al and Ti are studied in both the Altarelli-Feruglio A4 
model and the ^4 model of Ding. The rates of these lepton flavor violation process for 
the inverted hierarchy neutrino mass spectrum are enhanced considerably by the next to 
leading order contributions. For both models, we flnd that the /i — e conversion in Ti is 
within the precision of next generation experiments in all the parameter space considered, 
the fi — e conversion in Al should be observable at least in a very signiflcant part of the 
parameter space, whereas r 67, r fi^, r — ^ 3e and r ^ 3/x are below the expected 
future sensitivity. The detectability of ^ 67 and /x ^ 3e in near future depends on the 
models and the neutrino mass spectrum. We suggest that a comprehensive consideration 
of the lepton flavor violation processes is important to test and distinguish both discrete 
flavor symmetry models. 
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1. Introduction 

In the past decades, one of the most striking developments in particle physics beyond the 
standard model (SM) is the experimental establishment of neutrino masses and the large 
mixing property which is quite different from the small mixing in the quark sector. A 
global fit to the current neutrino oscillation data from the solar, atmospheric, reactor and 
long baseline neutrino experiments gives the following 3a interval for the mixing parameters 
[1,2] 

sin^ = 0.304t°;°^^, sin^ 9^3 = 0.50^^;^^, sin^ 613 < 0.056 

Am^i = 7.65tofQ x 10"^ eV^ Am^^ = ±2.40tg;^^ x 10"^ eV^ (1.1) 

For the overall scale of the neutrino and the mixing angle 613 currently only upper limit 
exists. But considerable progress is expected from future double beta decay [3] and reactor 
neutrino oscillation experiments [4, 5]. It is remarkable that such a mixing pattern is 
very close to the so-called tri-bimaximal (TB) mixing [6], which leads to sin^ ^12 = 1/3, 
sin^023 = 1/2 and sin^ ^13 = 0. Thus it is a important subject to realize the TB mixing 
from the theoretical viewpoint. 

Recently it is found that the TB mixing appears naturally in models with discrete 
flavor symmetry such as A4, T' and ^4, and many models based on these flavor symmetries 
have been proposed [7, 8, 9, 10, 11, 12, 13, 14]. We note that the T' model usually has the 
same structure as that of the ^4 model in the lepton sector. In these models, the flavor 
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symmetry is realized at a higher energy scale, the lepton fields transform nontrivially 
under the corresponding flavor group, and the flavor symmetry is spontaneously broken 
by a set of flavons along appropriate directions to provide a realistic description of the 
lepton masses and mixing angles. To be concrete, in this work we will concentrate on the 
Altarelli-Feruglio A4 model [9] and the S4 model of Ding [12], which represent a typical 
class of flavor models. The structure of the model is tightly constrained by the flavor 
symmetry, as a result only few parameters are involved at leading order (LO), and the 
models are rather predictive. The common features of both models are as follows: they 
are supersymmetry (SUSY) models and the neutrino masses are generated via the type-I 
Seesaw mechanism. At LO the charged lepton mass matrix is diagonal, the light neutrino 
mass matrix is controlled only by two complex parameters, and it is diagonalized by the 
TB mixing matrix exactly. Despite of the tight constraints on the model parameters, and 
in particular the fact that there is only one real parameter left after taking into account 
Am2i and Amg^^ measured from the neutrino oscillation experiments, it has been shown 
that the A4 models can reproduce the observed matter-anti matter asymmetry for natural 
values of the parameters in a quite satisfactory way [15, 16, 17, 18, 19]. 

Since this kind of A4 and ^4 models are so attractive and predictive, it is highly 
desirable to verify or exclude these model experimentally. There is no double that the 
precise measurement of 9is is a crucial test of the models. One of the implications of 
the observation of neutrino oscillation is the possibility of measurable branching ratio for 
charged lepton flavor violating (LFV) decays. While the LFV processes are still highly 
suppressed in the SM, predictions of the SUSY Seesaw for these rare decays are much more 
enhanced, as these processes are suppressed by the SUSY scale rather than the Planck scale. 
Furthermore, as different models obtain large neutrino mixing angles through different 
mechanisms, their predictions for the LFV decays can be very distinct. In the past years, 
LFV processes have been explored extensively from experiments, the upper bound of the 
branch ratio would be improved considerably in future, which could strongly constrain the 
new physics beyond the SM. Consequently, LFV rare decays may provide a way to test the 
A4 and 5^4 flavor symmetries. 

In this work, we shall investigate the predictions for the various LFV processes in 
the Altarelli-Feruglio model and the Ding's ^4 model, assuming the so-called minimal 
supergravity (mSUGRA) boundary conditions. Although our analysis is performed for the 
two specific models, it has generic features which are universal to the models based on the 
discrete flavor symmetry. We find that these models are so predictive that the branching 
ratios of the LFV processes are closely related to light neutrino mass, and there are some 
characteristic relations between various LFV branching ratios. Present work is different 
from the previous studies of the LFV decay branching ratio within the SUSY Seesaw 
framework. Since the low energy data allow to reconstruct only partially the high energy 
couplings of the theory, they generally inversed the Seesaw formula following the procedure 
of Casas and Ibarra [20], then they carried out the Monte Carlo studies by scanning the 
unknown right handed neutrino mass spectrum and the angles and phases of the inversion 
matrix in order to present scatter plots of the rare branching ratios [21, 22, 23]. Another 
method is the top-down approach, where the neutrino Yukawa coupling matrix and the 
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right handed neutrino mass matrix is determined by a specific SUSY grand unified theory 
[24]. 

The paper is organized as fohows: in the second section we briefly review the LFV 
within the framework of the SUSY Seesaw. Then we give a concise introduction to the 
Altarelli-Feruglio model and the Ding's S4 model in section 3. Our analysis of LFV decay 
branching ratios for these two interesting models is presented in section 4 and section 5 
respectively. Finally we summarize our results and draw the conclusions. 



2. LFV in SUSY Seesaw 

SUSY is a well motivated possibility for new physics [25, 26], and the supersymmetric See- 
saw mechanism can accommodate simultaneously tiny neutrino masses and large hierar- 
chy between the electroweak scale and the high Seesaw scale without serious fine-tunings. 
In this framework the particle content of the minimal supersymmetric standard model 
(MSSM) is extended by including three gauge singlet right handed neutrino superfields 
(z=l, 2, 3). Imposing R— parity conservation, the leptonic part of the superpotential is 
thus given by 

Wiep = elYe.jLjHd + u^Y^ijLjHu + ^ ^y!M^^j^y■ (2.1) 

where Li denotes the left handed lepton doublet, is the right handed charged lepton, 
Hu and Hd are the hypercharge +1/2 and —1/2 Higgs doublets respectively. Yg and 
Yjy are the charged lepton and neutrino Yukawa coupling matrices respectively, and M^^ 
represents the 3x3 major ana mass matrix for the heavy right handed neutrinos. In general, 
the overall scale of M^^, which we shall denote by M^^j, is much larger than the electroweak 
scale. At energies below M^^j, the theory will be well described by the following effective 
superpotential 

W[J/ - e|Ye^JL,Ha + ^A^..,(L,i/„)(L,if„) (2.2) 

where the light neutrino effective mass matrix is given by the well-known Seesaw 
formula 

M,^-YlM-'Y,vl (2.3) 

where Vu is the vacuum expectation value of the neutral component of the Higgs doublet 
Hu. In nature SUSY must be a broken symmetry, we thus introduce the soft SUSY breaking 
(SSB) Lagrangian Cgofti which could be new sources of flavor violation. The leptonic part 
of Csoft has the following form, 

-Csoft = {ml)ijLlLj + {io4)ije'iej + (m|)^^-?*?^- + ]^{Ae)ije'iLjHd 

+ {A,),,u:L,Hu + h.c] (2.4) 



where L^, and Ui are the supersymmetric partners of the left handed lepton doublet, right 
handed charged lepton and right handed neutrinos respectively. m|, m| and m~ are the 
corresponding soft mass matrices squared, Ag and Ajy are the charged lepton and neutrino 
soft trilinear couplings respectively. In general, above SSB terms can have arbitrary flavor 
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structures which would induce unacceptably large flavor violating effects. The simplest 
solution to this SUSY flavor problem is to assume a flavor-blind SUSY breaking mediation 
mechanism, which will generate flavor universal SSB terms at some high scale. In the 
present work, we will restrict ourselves to the so-called minimal supergravity scenario 
(mSUGRA). It assumes that, at the GUT scale, the slepton mass matrices are diagonal 
and universal in flavor, and that the trilinear couplings are proportional to the Yukawa 
couplings 

where mg is the common scalar mass, and Aq is the common trilinear parameter. In 
addition, there are still three parameters characterizing the mSUGRA: the common gaugino 
mass Ml/2, tan/3 = Vu/vd and the sign of the Higgs mixing parameters fi, where Vd is the 
vacuum expectation value of the neutral component of the Higgs doublet H^- However, 
these SSB terms will not be universal at the weak scale. Since the Yukawa coupling matrices 
Ye and Y^, can not be simultaneously diagonalized, non-vanishing off-diagonal elements 
would be generated in the SSB mass-squared matrices and the trilinear couplings due to 
renormalization effect. The one-loop renormalization group equations for m|, m| and Ag 
are as follows [22, 27], 



+(YtY, + YtYe)m| + 2At A, + 2AtAe - Qgl^^ - ^glmi^ - ^gfS 



dmi 1 



(2m| + 4m^jYeYt + 4Yem|Y| + 2YeYtm| + 4AeAt 



= ^2{^em^^dy\ + YeYt) + 5YtYe + yJ Y. - 3^1 - lal] 

+Ye[Tr(6AdYj + 2AeYt) + 4YtAe + 2Yt A, + G^i^a + y 5?Mi] } (2.6) 

where fi is the renormalization point, and we have defined 

S = ml^- ml^ + T¥[m| - m| - 2m| + m|+ m|] (2.7) 

Here m| cind m?- are respectively the SSB mass-squared matrices for the supersym- 
metric partner of the left handed quark doublet, right handed up type quark and right 
handed down type quark. In the phenomenological studies it is convenient to work in the 
leptonic basis where both the charged lepton mass matrix and the mass matrix of the right 
handed neutrino are real and diagonal. In the leading-logarithmic approximation with uni- 
versal boundary conditions Eq.(2.5), the off-diagonal elements of the above SSB slepton 
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mass matrices and trilinear coupling are given by [28, 29] 

(m|),^, ~ (2.8) 

where the hat denotes the basis in which the right handed neutrino and the charged lepton 
mass matrix are real and diagonal, and the factor L is defined as 

Lii = log(^)% (2.9) 

Here Mi is the ith heavy right handed neutrino mass, and Mq is the GUT scale typically 
equal to 2 x 10^^ GeV. Eq.(2.8) shows that, within the type-I Seesaw mechanism the off- 
diagonal elements of the right handed charged slepton matrix approximately don't run in 
the leading-log approximation, and the running of the trilinear parameters Ag is suppressed 
by the charged lepton masses. At low energy, the flavor off-diagonal correction (mj)ij 
induces the LFV processes such as k Ij^, k 3lj and k — Ij conversion in nuclei, with 
i > j = 1, 2, 3. A very useful tool to treat analytically the complicated expression for the 
branching ratio is the mass insertion approximation, where the small off-diagonal elements 
of the soft mass matrix are treated as insertions in the sfermion propagators in the loop. 
Then we obtain the compact expression for the branching ratio of the charged lepton LFV 
radiative decay 

Br{k ^ Ijl) ^ Br{k ^ /,^>,^^)^^^|(m|),,f tan^ /? (2.10) 



G],mi 

where Gi? is the Fermi coupling constant and is the characteristic mass scale of the 
SUSY particle in the loop. It has been show that the expression (2.10) with (m~)^j given 
by Eq.(2.8) represents an excellent approximation to the exact renormalization group result 
if one sets [27] 

ml {).^mlMli^{ml + Q.QMl/^f (2.11) 

The LFV processes U ^ Ij^, k 31 j and k — Ij conversion in nuclei occur via 7—, Z— 
and Higgs-penguins as well as squark/slepton box diagrams. It has been shown that the 
contribution from the 7— penguin diagrams are almost undistinguishable from the total 
rates in the universal mSUGRA scenario, although the contribution of Higgs-penguins 
becomes large in the region of large tan /? ^ 60 and light Higgs boson mass^ 100 GeV, 
which is not allowed by the current experimental lower bounds on the MSSM particle masses 
[34]. In this approximation, the branching ratio for the trilepton decays is approximately 
given by [21, 30] 

Br{k -> 3/,) ^ ^(log ^ - ^)Br{k - Ijl) (2.12) 
37r mf 4 

For fji — e conversion, the 7— penguins dominance gives 

CR(,N ^ eN) = ^ ZZt„\F(erBHj. ^ e-,) (2.13) 

-L cap iZTT i cap 
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BR(r 




4.5 X 10~^ 


10 ^ 


BR(r ^ 
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X in~^ 


-L VJ 


BR(/i ^ 


eee) 


1.0 X 10-12 


10-14 


BR(r ^ 




3.2 X 10-^ 


10-9 


BR(r ^ 


eee) 


3.6 X 10-^ 


10-9 


CR(/iTi 


^eTi) 


4.3 X 10-12 


10-18 


CR(/^A1 


^eAl) 




10-16 



Table 1: Present and upcoming experimental limits on various lepton flavor violation processes 
[34]. 

where Tcap is the measured total muon capture rate, Z is the proton number of the nucleus, 
Zf^ff is the effective atomic charge obtained by averaging the muon wavefunction over the 
nuclear density, and F{q^) denotes the nuclear form factor at momentum transfer q [31]. 
In this work, we will consider /i — e conversion in two nuclei ^^Al and 22^1. 13AI would 
be used by the proposed Mu2e experiment at Fermilab [32], Z^ff = 11.5, the values of 
the relevant parameters are Fi^cf' 2^ ~'^^) — 0.64 and Pcap — 

4.64079 X 10-19 GeV. For 

the nucleus "^Ti used by the proposed PRIME experiment at J-PARC [33], Z^ff = 17.6 , 
F{q^ -ml) 0.54 and Tcap = 1.70422 x 10-^^ GeV. 

Experimental discovery of LFV is one of smoking gun signatures of new physics beyond 
the SM, thus many experiments have been developed to detect the rare LFV processes. 
The present and projected bounds on these processes are summarized in Table 1. In the 
following, we shall investigate in details the predictions for the rare processes li Ij^, 
k 31 j and k — Ij conversion in ^3 Al and 22"^^ model and the S4 model of Ding, 

and discuss the possibility of testing these discrete flavor symmetry models by LFV. 

3. Flavor models with A4 and S4 symmetries 

In this section, we will briefly present the Altarelli-Feruglio A4 model [9] and the ^4 model 
of Ding [12], the next to leading order (NLO) corrections to the lepton mixing parameters 
will be discussed in detail, especially the allowed region of the reactor neutrino angle 9is 
and the Jarlskog invariant J in the lepton sector [35]. 

3.1 AF A4 model 

The AF model is a typical supersymmetric flavor model with A4 symmetry, and the aux- 
iliary symmetries are Zs and the Froggatt-Nielsen symmetry U{1)fn- The U{1)fn is to 
reproduce the observed charged lepton mass hierarchies, and the Zs symmetry is to guar- 
antee that the neutrino and charged lepton couple with different flavons at LO, so that the 
A4 symmetry is broken down to Z3 and Z2 subgroups in the charged lepton and neutrino 
sectors respectively at LO, this misalignment of symmetry breaking is exactly the origin of 
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TB mixing. A4 is the group of the even permutation of four objects, geometrically it is the 
invariant group of the regular tetrahedron in 3-dimensional space. It has three independent 
one-dimensional representations 1, V and and one three-dimensional representations 3. 
The multiplication rule for two triplet representations is the following 

3 3 = 1 e 1' e 1'' e 3^ e 3a (3.1) 

If we denote the two triplets as (ai, a2, as) and (61, 62, ^3)5 the singlets and triplets contained 
in the product are given by 

1 ^ aibi + a2bs + asb2 
1^ ^ a^bs + aib2 + a26i 
1^' ~ a262 + ^163 + asbi 

3s ^ (2ai6i - a2b3 - asb2, 2asbs - aib2 - a2bi,2a2b2 - aibs - a^bi) 

3a ^ (^263 - ^362, ^162 - ^261, asbi - aibs) (3.2) 

The fields in the model and their transformation properties under the flavor group are 
listed in Table 2. At LO the flavon fields develop the following vacuum expectation values 
(VEVs), 

(^t) = (^t,0,0), {e) = Ve, 

M-iVs^Vs^Vs) (3.3) 
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Table 2: The fields of the A4X Z^x U{1)fn model and their representations, where uo is the third 
root of unity co = e*^^/^. 

At LO the neutrino masses are generated by the following superpotential 

= y:.{jy''£)hu + {xa^ + XAO(i^V^) + ^^(zyVV^) (3.4) 

Here and in the following, we denote an invariant under A4 by a parenthesis (...). After 
the electroweak and flavor symmetry breaking, we have 

Wr, = yj,{u^jyi+u^jy3-\-iy3Jy2)vu+XAV^{uf+2u^jy^)+2x 

(3.5) 
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where Vu = {hu) is the VEV of the up type Higgs, and one can always set yi, to be real 
by a global phase transformation of the lepton doublet field. The first term in Eq.(3.5) 
contributes to the neutrino Dirac mass matrix 



/l o\ 
1 

[oioj 



(3.6) 



The remaining terms lead to neutrino Majorana mass matrix 



/ A+^B -\B 



Mm 



\ 



3 

A- 



-IB 
A-\B 

IB 



\ 



(3.7) 



where A = 2x^V^ and B = QxbVs- The right handed neutrino mass matrix Mm can be 
diagonalized by a unitary matrix U as usual 



U^MmU 



diagi\A + B\,\A\,\~A + B\) 



The unitary matrix U is given by 

U = UtbU^ 

where Utb exactly is the well-known TB mixing matrix 

/ ,/2 ^ Q \ 

Utb = 



(3.8) 



(3.9) 



2 J_ 
3 

111 

V6 V3 V2 



(3.10) 



L J L / 

\ V6 V3 V2 / 



and = diag(e-^^i/^e-^^2/2^g-i(/>3/2) jg ^ matrix of phase with = arg(A + B), 02 = 
arg(74) and 03 = arg(— A + B). The hght neutrino mass matrix is given by the type-I 
Seesaw formula 

= -MIM]^^Md (3.11) 
It is diagonahzed by the unitary transformation Uy^ i.e., 



rrijyUjy diag(mi, 7712, ^3) 
where Uiy = iUTsU^, and mi^2,3 are the hght neutrino masses 



mi = 



\A + B\ 



m2 = 



2 2 
\A\ ' 



ms = 



(3.12) 



(3.13) 



It is obvious that the AF model is strongly constrained, the neutrino part of the Lagrangian 
depends on only three parameters: y^^, A and the latter two parameters are in general 
complex numbers. Concerning the light neutrino at low energy, only two parameters A/y^ 
and B /y1 are involved. Imposing the constraints of Am2i and Amg^^ measured from the 
neutrino oscillation experiments, only one real parameter is left, and it is usually chose to 
be the lightest neutrino mass. Detailed studies showed that the flavon VEVs Vr, and 
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Vs should be approximately of the same order (9(A^A) [9, 16], where Ac is the Cabibbo 
angle, and A is the cutoff scale of the theory. The neutrino mass spectrum can be either 
normal hierarchy (NH) or inverted hierarchy (IH), and the lightest neutrino mass is tightly 
constrained as follows [16, 17] 



0.0044 eV < mi < 0.0060 eV, NH 
O.OlTeV < ma, IH 



(3.14) 



The above LO predictions for the lepton masses and mixing parameters could be cor- 
rected by both the higher dimensional operators in the superpotential Wi, compatible with 
the symmetry of the model and the shift of the vacuum alignment. Including the NLO 
operators, the charged lepton masses are described by the following superpotential. 



+yeij^e''{£(pT^T)0'^hd 



y^i^^i\t^T^T)'Ohd + yri^T\l^T^T)"hd (3.15) 

where 5lpt denotes the shifted vacuum of the flavon lpt- The first line represents the LO 
contributions, which leads to a diagonal charged lepton mass matrix, and the latter two 
lines are the NLO corrections. Since {{(^t^t)^s) ~ 0' 0) which has the same alignment 
directions as {(fr)^ consequently the contributions of ^gi, V/j,! and yri terms can be absorbed 
into the LO results. Taking into account {S(Pt) — (1,1, 1)5Vt [9], the charged lepton mass 
matrix at NLO is given by 



mi 



^6^2 A ye ^2 A i^eA2 A ^ 
„ y^m^ „ Ve_Vr Ve SVt 

\ yr A A y^ A / 



Vd 



(3.16) 



where we have redefined Vr + SVt ~^ Vt- The charged lepton mass matrix m^ can be 
diagonalized by the bi-unitary transformation, i.e., ulcm(Ui — diag(me, m^, m^), and the 
unitary matrix Ui approximately is 

/ 



1 

SVt 
' Vt 

5Vt 
' Vt 



/ 5Vt \ ^ 
V Vt ^ 

1 

5Vt 
' Vt 



/ SVt ' 

V Vt > 
/ SVt ' 

V Vt > 

1 



(3.17) 



(3.18) 



The NLO corrections to the neutrino Dirac mass terms are as follows, 

where (•••)3s(a) denotes the 'is(A) product of the two triplet representations, which can be 
read directly from the multiplication rule Eq.(3.2). Substituting the VEVs in Eq.(3.3), we 
obtain the NLO correction 5Md 

hyA \ 

5Md= -VA + VB 
-VA -VB 



V 



Vt 



(3.19) 
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The Major ana mass matrix of the right handed neutrino are corrected by the following 
terms at NLO 



(3.20) 



where S^, 6^ and 5(ps denote the shifted vacuum of ^, ^ and (ps respectively. Absorbing 
the above corrections partly into the LO results, then SMm can be parameterized by four 
independent parameters as follows 



SMm = 



^ 2xd xb - xe xc ^ 

xb - Xe Xc -xd 
y Xc -Xd Xb + 2x£; j 



Vt 

T 



B 



(3.21) 



Therefore at NLO the light neutrino mass matrix is given by mjy = —{M'£ + 5M'£){Mm + 
SMm)~^{Md + SMd)- Performing the standard perturbation diagonalization, niu can be 
diagonalized by the unitary transformation Ul as U'J-nii^Ul = diag(m^, m3). To first 
order in the expansion parameter Vr/A = e, the light neutrino masses wi'^^ 2,3 



1712 = 



A + B 2{A + Bf 



[AvaA + AvaB + {xb + XC- 2xD - 2xE)y^B]e 



A 



+ 



A^ 

yvvl 



{xb + xc)e 

-,[-AyAA + AvaB + {xb + xc + 2xd + 2xE)yuB]e (3.22) 



2 _ 1 * 

3 V3^12 



73 + 



^13 + ^^23 ^ 



A-B 2{A-BY 
The mixing matrix Ul is given by 
/ 

\~7g~ 73^12 + 73*13 73 " 76*12 + ^S23 - + 75^23 / 

The small parameters S12, S13 and S23 are 
1 



"xJe v)3*12 V2*13 ^ ^«12 ^^23 ^ yg^lS + ^3*23 



(3.23) 



Sl2 = 



■[A*[AyAA + 2yBB + {-2xd + XE)y,B]e + {A + B)[4y*^A* 



V2{\B\'^ + AB* + A*B)y^ 
+ 2y%B* + {-2x*D + x*E)y,B*]e*} 

si3 = 4^(^AB*\A*B)yS ~ ~ ^ ^^^"^ ~ ^""^^'^^ ^ ^ ^^^^^^ ^ ^^"^^ ~ 

S23 = ^^AB. + A*B-myS ^^^'^ " + ^^^^'^^^ ~ + 35^y.)e*} (3.24) 
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Taking into account the NLO corrections to the charged lepton mass matrix, the leptonic 
PMNS matrix becomes UpMNS — U^U^, consequently the parameters of the lepton mixing 
matrix are modified as 



e3\ 



1 



+ 



sin^ ^21 



1 

3 ~ 3 



3\/2(AB* + A*B - |B|2)y^ 
2 



I _ B{A* - + 3(XB - xc)yu\e + + B)\\yl + 3(x^ 

\b{A* - B*){2yB + 3xEy.)e - B*A{2yl + Zx*Ey,)e*] 
U2A* + B*)[AyAA + 2yBB + {-2xd 



'i{\B\'^ + AB* + A*B)y^ I 



+XE)yuB]e + {2 A + B)[4.y\A* + 2y|jB* + (-2x|, + x|;)y,S*]£*} 
1 1,5Vt ,5Vt... 1^12 



•2/1 ± ,u t'x /" ^^T N*-! 

sm ^23 = 2 + 2 + (-^) ] - 12(^5* + A*B)y, 



{ [4yB + 3{XB - xc)y^]e + [4y^ + S{x% - 



-xc)yu]e*^ - 



\B 


2 


6{AB* + A*B - 


B 





(3.25) 



It is obvious that all the three mixing angle receive corrections of the same order of mag- 
nitude about A^, which is allowed by the current neutrino oscillation data in Eq.(2.1). In 
particular, the reactor mixing angle ^13 becomes of order after including the NLO cor- 
rections. In order to see clearly the behavior of the mixing parameters after including the 
subleading contributions, we perform a numerical analysis. The LO and NLO parameters 
Ha^ Vb^ xb^ xc^ xd and xe are treated as random complex number with absolute value 
between and 2, A and B in the right handed neutrino mass matrix Eq.(3.7) are random 
complex number with absolute value in the range of 10^^ — 10^^ GeV, and the parameters 
5Vt/Vt and V^/A have been fixed at the indicative value of 0.04. Furthermore, we require 
the oscillation parameters to lie in their allowed 3cr range given in Eq.(2.1). The results are 
presented in Fig. 1, where the Jarlskog invariant J of CP violation in neutrino oscillation 
is defined as 



J Im[{UpMNs)el{UpMNs)i22{UpMNs)l2{UpMNS^ 



(3.26) 



As we can see from this figure, in AF model 6is should be much smaller than the present 
upper bound, and it prefer to lie in the la range. Our analytical estimates are confirmed. 



3.2 5^4 model of Ding 

The total flavor symmetry of this model is S4 x Zs x Z4 [12], where the auxiliary symme- 
try Z3 X Z4 is crucial to eliminate unwanted couplings and to insure the needed vacuum 
alignment. It is remarkable that the realistic pattern of fermion masses and flavor mixing 
in both the lepton and quark sector have been reproduced in this model, and the mass 
hierarchies are controlled by the spontaneous breaking of the flavor symmetry instead of 
the FN mechanism. The matter fields and the flavons of the model and their classification 
under the flavor symmetry are shown in Table 3, where the quark fields have been omitted. 
In this section, the convention for the group theory of ^4 is the same as that in Ref.[12]. 
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Figure 1: Scatter plot of sin O23 and the Jarlskog invariant J against sin Ois for both normal 
hierarchy and inverted hierarchy neutrino mass spectrum in AF model. 
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Table 3: The transformation rules of the matter fields and the flavons under the symmetry groups 
6'4, Z3 and Z4 in the ^4 model of Ref. [12]. uo is the third root of unity, i.e. uo = = ( — l+i\/3)/2. 



Explicit calculation demonstrated that the LO vacuum alignment is as follows [12], 



(^) = (o,y^,o), (x) = (o,y^,o) {0 = Vc 



(3.27) 



In this model, the charged lepton masses are described by the following superpotential 

+ ^ e'^((i'x)3i(w)3i)ii/id+ |^/^'^(%x)3.)i2/id+ + ... (3.28) 

Taking into account the vacuum alignment in Eq.(3.27), we find that the charged lepton 
mass matrix is diagonal at LO, and the charged lepton masses are given by 



me 



m. 



rrin- 



(3.29) 
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where is the hnear combination of y^i {i — 1 — 10). At LO the superpotential contributing 
to the neutrino mass is as fohows 



(3.30) 



The neutrino Dirac and Majorana mass matrices can be straightforwardly read out as 



/ 26 a-h a-h\ 

a — b a + 2b —b 
ya-6 -b a + 2b J 



M \ 
Vu, Mm = I M 
M y 



(3.31) 



where we denote a — yj^iI^/A and b — yu2V(j,/A, both of them should be of order A^. The 
heavy right handed neutrino mass matrix Mm can be diagonalized as follows 



U^MmU = diag(|M|, |M|, \M\ 



(3.32) 



It is obvious that the right handed neutrinos are exactly degenerate, and the unitary matrix 
U is 

/I \ 



U 



^ ^ 



(3.33) 



Here Ua — diag(e~'^'^/^, e""^*^/^, ze~^^/^), and a — arg(M) is the complex phase of M. The 
hght neutrino mass matrix is given by the Seesaw relation m^^ = —M^M^M^^ which is 
exactly diagonalized by the TB mixing matrix 

U^myUy diag(mi,m2,m3) 



g -^2/2 ^-ioLz/2\ 



The unitary matrix Uy is written as 

C/, = [/TBdiag(e-^"i/2^ 
where the phase ai,2,3 and the light neutrino masses mi^2,3 are given by 

ai = arg(-(a - ^bf/M), a2 = arg(-4a^/M), as = arg((a + ^bf/M) 



(3.34) 



(3.35) 



mi 



(a - 36)^ 



m2 — 4|a 



ms = \{a + ^bf 



K 
\M\ 



(3.36) 



Similar to the AF model, this model is very predictive, there are only three independent 
parameter a, b and M in the neutrino sector. The neutrino mass spectrum can be normal 
hierarchy or inverted hierarchy as well. Taking into account the mass square difference 
Am^i and Aml^ measured in the neutrino oscillation experiments, we obtain the following 
limit for the lightest neutrino mass 



mi > 0.011 eV, NH 
ms > 0.0 eV, IH 



(3.37) 



In the following, we briefly discuss the NLO corrections, please read Ref.[12] for details. 
After including the higher dimensional operators allowed by the symmetry at NLO, the 
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off-diagonal entries of the charged lepton mass matrix become non-zero and of the order 
of the diagonal term in each row multiplied by 6, which parameterizes the ratio VEV/ K 
with order O(A^). Therefore the charged lepton mass matrix can generally be written as 



mi 



^3^ 



21^ 



23^ 



m 



31^ 



m 



32^ 



(3.38) 



33 / 



where the coefficients mlj{i,j =1, 2, 3) are order one unspecified constants. The matrix 
m\mi is diagonalized by setting I to U^l^ where Ui approximately is given by 



/ 1 AsY Aey\ 



''22 

1 



(3.39) 



^31^ 

^33 ^33 

The NLO correction to the Majorana masses of the right handed neutrino arises at order 
1/A, the corresponding higher dimensional operator is / whose contribution 

can be completely absorbed into the redefinition of the mass parameter M. Consequently 
the right handed neutrinos are highly degenerate, and we only need to consider the NLO 
corrections to the neutrino Dirac couplings as follows 



uhu + ^{ty''iv)i2(hu + ^(i^" 



(3.40) 



A2 -"^^^-^ ' A2 

A part of the above corrections can be absorbed into the LO results, then the NLO cor- 
rections to Md can be parameterized as 



5M 



D 







be ae — be 



—be ae 
yae + be —be 



be 




(3.41) 



where the magnitudes of d and b are expected to be of the same order as those of a and b. 
In first order of e, the parameters of the lepton mixing matrix are modified as 

1 1 



Ue3\ = 



V2 6(|a|2 + 9|6|2)(a6* + a*6) 

^ ^*_ ^ 



[(a -t- 36)^(a*a* + 66*6*)e* - (a* - 36*)^(aa + 6bb)e\ 



sm U12 
.2. 



sin 023 



22 
1 
3 
1 



33 



m: 



21 



m: 



22 



777,33 ^^22 



+ 



2 2(|a|2 + 9|6|2)(a6* + a*6) 



a6(a*a* + 66*6*)e* + a%%aa + 6bb)e 



+ 



1 [m^ 



32 



L777: 



33 



j^ ^32 . 



777: 



33 



(3.42) 



As we can see, NLO contributions introduce corrections of order to all mixing angles, 
and the reactor angle 9is becomes non-zero. Then we turn to a numerical analysis. Since 
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a, 6, a and h are expected to be of order A^, they are treated as complex numbers with 
absolute value between 0.01 and 0.1, the parameters m^^^jm^^^ ^31/^33 ^^id ^32/^33 in 
the charged lepton mixing matrix Un are taken to be complex numbers with absolute value 
in the range of 0-2, the heavy neutrino mass \M\ is allowed to vary from 10^^ GeV to 10^^ 
GeV, and the expansion parameter e is set to the indicative value 0.04. The correlation for 
sin^ 023 — sin^ ^13 and J — sin^ ^13 is showed in Fig. 2, It is obvious that rather small ^13 is 
favored, which is consistent with our theoretical analysis. 




Figure 2: Scatter plot of sin^ d^z and the Jarlskog invariant J against sin^ ^13 for both normal 
hierarchy and inverted hierarchy neutrino mass spectrum in the ^'4 model of Ref.[12]. 



4. Predictions for LFV in AF model 



We first study the branch ratios of the rare LFV processes at LO of the model. Since 
the number of independent parameters is rather small in this case, these branching ratios 
are closely related to the neutrino oscillation parameters. As has been shown in section 
3.1, The charged lepton mass matrix is diagonal at LO, and the right handed neutrino 
mass matrix Mm is diagonalized by the unitary matrix U — UtbUcJ), consequently in the 
base where both the right handed neutrino and charged lepton mass matrices are real and 
diagonal, the neutrino Yukawa coupling matrix is given by 



J_^-i(/)2/2 J_^-i(/>2/2 J_^-i(/>2/2 

Vs Vs Vs 

J_^-i(/)3/2 1 



v 







g .0i/2 ^ 

-i4>2/2 



V6 



(4.1) 



As a result we can straightforwardly obtain the off-diagonal elements of Y^LYi^, which is 
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directly related to the LFV branching ratios via Eqs.(2.8), (2.10), (2.12) and (2.13) 

(YtLY,)i2 = (YtLY,)2i = \yl\^ — 

6 mi 

(YtLY,)i3 = (YtLY,)3i = ^y^ln^ 

6 mi 



(Yt LY,)23 = (Yt LY,)32 = ^-yi In (4.2) 



1 2i mim\ 

-y^ In 3- 

mi 



Obviously the LFV processes are tightly related to the light neutrino mass, and the branch- 
ing ratio is proportional to y^. We notice that the branching ratio is independent of the 
grand unification scale Mq^ and it tends to zero if the neutrino mass spectrum is degen- 
erate. Whereas the branching ratio could become considerable large if the neutrino mass 
spectrum is strongly normal hierarchy or strongly inverted hierarchy. It is remarkable that 
we have (YjLY^y)i2 = (YjLY^y)i3, which is related to the fi — r symmetry of the light 
neutrino mass matrix. As a result, the ratio of the branching ratios for r ^ 67 and ji 
approximately is 

^J^^ ^ ^ Br{r -> eUeT^r) ^ 17.84% (4.3) 

This ratio is almost independent of the SUSY breaking parameters. Given the present 
experimental bound Br{ii 67) < 1.2 x lO"-'^-'^ [34], this result implies that r — > 67 
has rates much below the present and expected future sensitivity [34]. We note that this 
ratio would be corrected drastically by the NLO contributions. For the IH neutrino mass 
spectrum, we have 



mi — y TTig + lAmg^l, 7712 — y + lAmg^^l + Am2]^, 

In^ = ll^/!!i±^!M^±^!!4x ^ ^^ii /4 4^ 

mi 2 ^ m§ + |Am|i| 2(m§ + |Am|i|) ^'^ 

As a result, in the case of IH spectrum, the LFV processes fi 67, r 67, fi — > 3e, 
T 3e and jj. — e conversion in nuclein generally have rather small branching ratios at LO, 
which should be smaller than those of the NH case. While for r — > /X7 and r — > 3//, the 
branching ratios for NH and IH spectrum should be approximately of the same order. 

The analytical results in Eq.(4.2) allow us to estimate the branching ratios of the LFV 
processes via the formula Eq.(2.10)-Eq.(2.13). For definiteness, we will present our results 
only for the well-known mSUGRA point SPS3 [36], which is taken as the reference example. 
The SPS3 point is in the co-annihilation region for the SUSY dark matter, and the values 
of the universal SSB parameters are as follows 

mo = 90 GeV, M1/2 = 400 GeV, Aq = GeV, tan /? = 10 (4.5) 

Since no suppression is expected for parameter that is unrelated to the breaking of the 
flavor symmetry, the coupling yjy should be of order 1. We set yi, to be equal to 0.5 in the 
following numerical analysis, the corresponding predictions for Br{£i ^jl)^ Br{£i 3£j), 
CR{ii — e,Al) and CR{ii — e,Tz) are shown in Fig. 3 and Fig. 4 for NH and IH spectrum 
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respectively. The uncertainties of both Am^i and Amg^^ are taken into account, they 
are allowed to vary within their 3a allowed range in Eq.(2.1). We can clearly see that 
the LFV processes r 67, r r ^ 3e and r ^ Sfi are below the present and 

future experimental precision at B factory for both the NH and IH spectrum. As has been 
stressed above, the predictions of Br{ii 67), Br{r 67), Br{ii 3e), Br{r 3e), 
CR{ii — e, Al) and CR{ii — e, Ti) for IH are really much smaller than those of NH case, the 
rates of r — > /i7 and r — > for NH and IH roughly have the same order. 



0.00440.00460.00480.00500.00520.00540.00560.00580,006C 



0.00440.00460.00480.00500.00520.00540.00560,00580.0C 



0.00440.00460.00480.00500.00520.00540.00560.00580.006C 



0.00440.00460.00480.00500.00520.00540.00560.00580.006C 



0.00440.00460.00480.00500.00520.00540.00560,00580,00 




0.00440.00460.00480.00500.00520.00540.00560.00580.006C 



Figure 3: Dependence of Br{li ijj), Br{£i 3£j), CR{ii - e,Al) and CR{ii - e,Ti) on the 
lightest neutrino mass mi in AF model for the normal hierarchy spectrum. The bands have been 
obtained by varying /S^m^i and Am^^ in their 3cr experimental range. The dashed and dotted lines 
represent the present and future experimental sensitivity respectively. There is still no upper bound 
for CR{fi — e, Al) so far. 





Figure 4: Dependence of Br{£i £jj), Br{£i 3£j), CR{fi - e,Al) and CR{fi - e,Ti) on the 
lightest neutrino mass ms in AF model for the inverted hierarchy spectrum. The bands have been 
obtained by varying Am2i and Am^i in their 3(j experimental range. The dashed and dotted lines 
represent the present and future experimental sensitivity respectively. 
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For the NH mass spectrum, numerical calculations demonstrate that the ratio between 
Br{r e^) and Br{fi 67) indeed is 17.84%, and Br{r fi'^) is about 2-4 times as large 
as Br {11 67). Therefore we conclude that the branching ratios oi ji ^ e7, r ^ 67 and 
T ^ ji^ roughly are of the same order, the same result has been obtained in Ref. [37]. It is 
notable that the current experiment limit Br{ii — > 67) < 1.2 x 10~^^ already constrains the 
model strongly, the lightest neutrino mass mi very close to the lower bound 0.0044 eV is 
disfavored. The rates of /x ^ 3e, fi — e conversion in Al and /i — e conversion in Ti are below 
the present upper bound, while they are above the expected future experimental sensitivity, 
consequently it is very promising to detect these three rare processes in near future. If this 
turns out to be true, the parameter space of the model would be strongly constrained. We 
plot the constraint on the parameters mi and imposed by the observation of /x ^ 67, 
/X ^ 3e and fi — e conversion in Al and Ti in Fig. 5. It is obvious that the branching ratios 
of these LFV processes are more sensitive to the parameter yy than mi. For y^ of order 1, 
we should be able to observe fi 67, fi ^ 3e and 11 — e conversion in Al and Ti at next 
generation experiments in the case of NH spectrum, in particularly the /i — e conversion 
processes, which is an even more robust test to the AF model. 




0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.5 1.0 1.5 
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Figure 5: The contour plots for 5r(/i 67), Br{fi 3e), CR{ii — e^Al) and CR{ii — e^Ti) in 
the parameter space of and mi for the NH spectrum in AF model, where Am2i and Amg^ are 
chose to be the best fit values 7.65 x 10~^ eV^ and 2.40 x 10~^ eV^ respectively. 



For the IH spectrum, the branching ratios for the jie and re involved LFV processes are 



indeed smaller than those of NH case, as has been stressed below Eq.(4.4). We notice that 
all the LFV processes are below both the present and future experiment sensitivity except 
the /J. — e conversion in Ti. However, the rates for fi 67, r — > 67, — ^ 3e, r 3e and 
II — e conversion in Al and Ti at LO are so small that they may be corrected considerably 
by the NLO terms. Therefore it is very necessary to include the NLO contributions and 
exploit the physical effects on the LFV observables. 

We perform a numerical analysis by treating all the NLO effective couplings yA^ Vb, 
xb, xc, xd and xe in Eq. (3.21) and Eq. (3.23) as random complex numbers with absolute 
value between and 2. The LO parameters A and B in the heavy right handed Majorana 
mass matrix Eq.(3.7) are taken to be random complex numbers with absolute value between 
10^^ GeV and 10^^ GeV, while and ^ has been fixed at the indicative values of 0.04. 
In order to compare with the leading order predictions, the coupling yy is set equal to 0.5 as 
well in the numerical analysis. The scatter plot of the LFV branching ratios vs the lightest 
neutrino mass for NH and IH spectrum are showed in Fig. 6 and Fig. 7 respectively. These 
plots display only the points corresponding to choices of the parameters reproducing Am2i, 
Am3]^ and the mixing angles within their allowed 3cr interval given in Eq.(2.1). 




Figure 6: Scatter plot of Br{U £^7), Br{£i 3£j), CR{fi - e,Al) and CR{fi - e,Ti) against 
the lightest neutrino mass mi in AF model for the normal hierarchy spectrum. The dashed and 
dotted lines denote the present and future experimental sensitivity respectively. There is no upper 
bound for CR^ji — e, Al) so far. 



We see that the LFV branching ratios for NH are modified slightly by the NLO contri- 
bution, the above discussions for NH case at LO still apply here, and the next generation 
experiments oi ji ^ 67, fi ^ 3e and ji — e conversion in Al and Ti are crucial tests to 
the AF model. However, for IH spectrum the rates of the processes ji 67, r 67, 
/X — > 3e, T — > 3e and /x — e conversion in Al and Ti are enhanced considerably. We note 
that Br{T — ^ ji^) and Br{T Sfi) are qualitatively the same as the LO results, this is 
because that r ^ jj.^ and r ^ Sfj. are related with the elements {Yt'LYjy)23, which are 
not suppressed at LO and are much larger than NLO corrections. It is notable that /i — e 
conversion in Ti is within the sensitivity of the next generation experiment for the whole 
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Figure 7: Scatter plot of Br{£i ijj), Br{ii 3ij), CR{ii - e,Al) and CR{ii - e,Ti) against 
the lightest neutrino mass ma in AF model for the inverted hierarchy spectrum. The dashed and 
dotted lines represent the present and future experimental sensitivity respectively. 



parameter space considered, the signals of /i ^ 67 and ji — e conversion in Al should be 
detected in near future in a very large part of the parameter space, and /i ^ 3e could 
be observed only in a marginal part of the parameter space. However, r 67, r fi^, 
T ^ 3e and r ^ 3/x are still below the expected future sensitivity. We note that the LO 
result ^^l^^g^^j — 17.84% is destroyed completely after including the NLO corrections, the 

variation of ^^^^^^^^j vs the lightest neutrino mass is presented in Fig. 8. 

In short summary, for both NH and IH spectrum of AF model, fi — e conversion in Ti 
can be observed in all the parameter space considered, /x — > 67 and fi — e conversion in Al 
should be observed at least in a very significant part of the parameter space, and fi ^ 3e 
may be observed on for NH spectrum. Of all the LFV processes, the jn — e conversion in 
Ti should be an even more robust one, since its sensitivity would be improved drastically 
in near future. 




6.002 0.004 0.006 0.008 0.010 'O-OO 0.05 0.10 0.15 0.20 



mi (eV) ms (eV) 

Figure 8: Scatter plots of Br{r ej)/Br{/j. ej) against the lightest neutrino mass for both 
NH and IH in AF model. 
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5. Predictions for LFV in S4 model 



Similar to section 4, we first study the predictions for the LFV branching ratios at LO 
of the model [12], then present the NLO corrections. From the discussion in section 3.2, 
we learn that the charged lepton mass matrix is diagonal, and the heavy right handed 
neutrinos are degenerate at LO . In the base where both the right handed neutrino and 
charged lepton mass matrices are diagonal and real, the neutrino Yukawa coupling matrix 
is = U^Yiy, where the unitary matrix U is given by Eq. (3.33). Straightforwardly we 
can calculate the hermitian matrix Y^LY^^ as follows 

/ 2|ap + 6|6p - 4|a||6| cos$ \a\^ - 3|6p + 2|a||6| cos$ |ap - 3|6p + 2|a 



yJly, 



V 



|ap-3|5p + 2|a||6|cos$ 2|ap + 6|6p + 2|a||6| cos$ \a\^ - 3\b\^ - 4\a\ 
lap - 3|5p + 2|a||6| cos$ |ap - 3|6|2 - 4|a||6| cos$ 2|ap + 6|6p + 2|a 

^ 2mi + 7712 ^2 — ^1 ^2 — ^1 ^ 

= 7712 — 7711 ^(7711 + 27712 + 37713) ^(7711 + 27712 — 377I3) 
y 7712 - 7711 ^(^1 + 27712 37713) |(771i + 27712 + 37713) J 

where $ is the relative phase between a and 6, and 7711^2,3 are the light neutrino masses 
given in Eq.(3.36). The off-diagonal elements of Y^LY^^ can be read out directly 



cos$ \ 
cos $ 
I cos $ J 



In 



\M\ Mg 
3^ ""Wl 



(YtLY,)i2 = (YtLY,)2i = (m2 - mi) 
(YtLY,)i3 = (YtLY,)3i = (m2 - mi) 
(YtLY,)23 = (YtLY,)32 = (mi + 2m2 



\M\ Mg 

^ W\ 
\M\ ^ Mg 



^ In — — 

3(mi + m2)vl \M\ 



Am^i 



3(mi + m2)vl 
Mg 



iMlln 



Mg 
\M\ 



(5.2) 



It is obvious that the off-diagonal elements of Y^^LY^^ are closely related to the light 
neutrino mass difference and the right handed neutrino mass |M|, and they increase with 
\M\. Interestingly we notice (YjLYiy)i2 = (YjLYiy)i3, which turns out to be true in 
AF model at LO as weh. The relation (yJlY^)i2 = (yJlY^) 13 seems universal in 
models which reproduce TB mixing at LO. Consequently we also have exactly Br{r 
e^)/Br{ji 67) 17.84% at LO, which is not satisfied anymore after including the 



subleading corrections, 
branching ratios oi ji - 



Since the factor 



is rather small for IH spectrum, the 



3(mi+m2) 

67, T 67, fi ^ 3e, r ^ 3e and fi — e conversion in Al and Ti 
should be small substantially at leading order. In the limit of degenerate light neutrino 
mass spectrum, the off-diagonal elements tend to be zero, the LFV processes would be 
highly suppressed. 

It is interesting to estimate the order of magnitude for the right handed neutrino mass 
\M\. Since the parameters a and b are expected to be of order A^, with this and using 
^/Am^ 2^ 0.05 eV as the light neutrino mass scale in the Seesaw formula, we obtain 



\M\ - 10^2 - lO'^GeV 



13r 



(5.3) 



We set \M\ to be equal to 1.0 x 10^^ GeV, the LO predictions for the branching ratios 
of LFV processes are plotted in Fig. 9 and Fig. 10 for NH and IH spectrum respectively. 



Mg 
\M\ 

(5.1) 
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where the SSB parameters are chose to be the SPSS point as well. We see that the LFV 
branching ratio is smaller the corresponding one of AF model, all the LFV processes are 
below the planned experiment sensitivity except — e conversion in Ti. Consequently the 
next generation experiment oi ii — e conversion in Ti is very important to test this 5^4 model, 
and the parameter space of the model could be tightly constrained by CR{/j. — e, Ti). The 
contour plot of CR{ii — e, Ti) in the mi — \M\ {rus — \M\) plane is showed in Fig. 11, where 
both NH and IH mass spectrum are considered. It is obvious that the contour curves move 
toward the left with the increase of the sensitivity. We notice that in this S4 model the 
signal oi fi — e conversion in Ti could be observed in a large part of the allowed parameter 
space. As we will demonstrate that in case of IH spectrum, CR{fi — e, Ti) receives relative 
large correction from the NLO contributions, as a result, the contour plot for IH should be 
taken with a grain of salt. 




Figure 9: Dependence of Br{£i £jj), Br{£i 3£j), CR{fi - e,Al) and CR{fi - e,Ti) on the 
lightest neutrino mass mi in the Ding 's S4 model for normal hierarchy spectrum. The bands have 
been obtained by varying Amli and Amg^ in their 3a allowed range. The dashed and dotted lines 
represent the present and future experimental sensitivity respectively. There is still no upper bound 
for CR{jii — e, Al) at present. 

Since the fie and re involved LFV processes are predicted to be suppressed at LO in 
the case of IH spectrum, the NLO contributions may be comparable to the leading ones. 
Therefore we should take into account the NLO contributions in order to reach a much 
more solid conclusion. In the following numerical analysis, the NLO effective parameters 
^2i/^225 ^31/^33 ^32/^33 Eq.(3.38) are treated as random complex numbers with 
absolute value between and 2. a, 6, d and b are taken to be complex random number with 
absolute value in the range of 0.01-0.1. The heavy neutrino mass |M| is fixed at 1.0 x 10^^ 
GeV in order to compare with LO results, and the expansion parameters e is set equal to 
0.04. The scatter plots of LFV branching ratios for the NH and IH spectrum are displayed 
in Fig. 12 and Fig. 13 respectively. 

We see that the LFV branching ratios can vary within a larger region than the LO 
results in Fig. 9 and Fig. 10, and the rates for /ie and re involved processes could be enhanced 
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Figure 10: Dependence of Br{£i £jj), Br{£i 3£j), CR{/j. - e,Al) and CR{/j. - e,Ti) on the 
lightest neutrino mass ms in the Ding's /S4 model for inverted hierarchy spectrum. The bands have 
been obtained by varying Am2i and Amg^ in their 3cr allowed range. The dashed and dotted lines 
represent the present and future experimental sensitivity respectively. 




Figure 11: The contour plot for CR{/i — e^ Ti) in the parameter space of the lightest neutrino mass 
and right handed neutrino mass |M| in the Ding's ^4 model, where Am2i and Am^^ are chose to 
be the best fit values 7.65 x 10~^ eV^ and ±2.40 x 10~^ eV^ respectively. 



by a factor a 10-100 for the IH case. However, Br{T fi^y) and Br{T 3/x) are only 
modified slightly by NLO contributions. It is remarkable that in this model the /i — e 
conversion in Ti is completely within the precision of the next generation experiment, the 
11 — e conversion in Al should be observable in a large part of the parameter space, the 
radiative decay /x ^ 67 can only be observed in a marginal part of the parameter space at 
near future experiment, and r 67, r fi^, fi ^ 3e, r ^ 3e and r ^ 3/x are below the 
expected future sensitivities. This fact can be used to distinguish this ^4 model from the 
AF model. 
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Figure 12: Scatter plot of Br{ii ijj), Br{£i 3ij), CR{ii - e, Al) and CR{ii - e, Ti) against 
the lightest neutrino mass mi in Ding's ^4 model for normal hierarchy spectrum. The dashed 
and dotted lines represent the present and future experimental sensitivity respectively. There is no 
upper bound for CR{fi — e, Al) so far. 
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Figure 13: Scatter plot of Br{£i £jj), Br{£i 3£j), CR{fi - e, Al) and CR{/j. - e, Ti) against 
the lightest neutrino mass ms in Ding's ^4 model for inverted hierarchy spectrum. The dashed and 
dotted lines represent the present and future experimental sensitivity respectively. 

6. Conclusions and discussions 

Recently some models with discrete flavor symmetry such as A4 and ^4 have been showed 
to be able to naturally produce TB mixing at leading order. Since only few parameters are 
involved at LO in these models, they are rather predictive. In this work, we have studied 
the LFV fi 67, r 67, r ^7, ^ 3e, r ^ 3e, r ^ 3/x and ji — e conversion in 
Al and Ti in the AF ^4 model [9] and the 6^4 model of Ding [12] within the framework of 
mSUGRA. 

At LO the branching ratio for LFV process is closely related to the light neutrino mass. 
For inverted hierarchy neutrino mass spectrum, fi 67, r 67, ^ 3e, r ^ 3e and fi — e 
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Table 4: Summary of predictions for LFV processes in AF model [9] and Ding's /S4 model [12]. 
The symbol ^J denotes the rare process is above the sensitivity of next generation experiments for 
the whole parameter space considered, / represents the LFV process should be observed in a very 
large part of parameter space in near future, ? for the process only observable in a marginal part of 
parameter space, and x denotes the process is below the sensitivity of next generation experiments. 



conversion in Al and Ti are highly suppressed at LO. After taking into account the NLO 
contributions, the LFV branching ratios for NH are corrected properly, and the allowed 
region becomes larger. Whereas the branching ratios of the fie and re involved LFV are 
enhanced drastically in the case of IH spectrum, and Br{T — > ji^) and Br{T are 
approximately the same as the LO predictions. 

Our predictions for the rare LFV processes in the two models are summarized in 
Table 4. From detailed numerical analysis, we find that for NH spectrum of AF model, 

67, fj. 3e and jj. — e conversion in Al and Ti are within the reach of next generation 
experiments (Fig. 6). While for IH spectrum fi — e conversion in Ti is above the expected 
future sensitivity in all the parameter space considered, fi ^ and fi — e conversion in 
Al could be observed by near future experiments in a very significant proportion of the 
parameter space, nevertheless the signal of 3e could be only detected in a marginal 
part of the parameter space (Fig. 7). Then we turn to the Ding's 5*4 model, for both 
NH and IH spectrum, we find that only fi — e conversion in Ti is within the precision of 
the future experiment, /i — e conversion in Al should be observed in a large part of the 
parameter space, while ji observable in a marginal part of the parameter space (Fig. 

12 and Fig. 13). In both AF model and the ^4 model of Ding, r ^ 67, r ^ /X7, r ^ 3e 
and T ^ 3/jL are below the sensitivity of next generation experiments. We conclude that 
ji — e conversion in Ti is a particularly robust test to the AF model and the Ding's ^4 
model. If it is really observed in future, it would be a great support to these discrete flavor 
models. 

We note that our consideration are not at all restricted to the AF model and Ding 's 
5^4 model, but could apply to a much wider class of theories. Models with discrete flavor 
symmetry can be strongly constrained or be excluded by existing or future LFV bounds, 
a combined analysis of LFV processes provide us a way to distinguish different models. 
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Note Added: near the completion of this work, papers that address similar issues appear 
[38, 39]. 
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